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Using tlie momentum-dependent MDI effective interaction for nucleons, we have studied the 
transition density and pressure at the boundary between the inner crust and liquid core of hot 
neutron stars. We find that their values are larger in neutrino-trapped neutron stars than in neutrino- 
free neutron stars. Furthermore, both are found to decrease with increasing temperature of a neutron 
star as well as increasing slope parameter of the nuclear symmetry energy, except that the transition 
pressure in neutrino-trapped neutron stars for the case of small symmetry energy slope parameter 
first increases and then decreases with increasing temperature. We have also studied the effect of the 
nuclear symmetry energy on the critical temperature above which the inner crust in a hot neutron 
star disappears and found that with increasing value of the symmetry energy slope parameter, the 
critical temperature decreases slightly in neutrino-trapped neutron stars but first decreases and then 
increases in neutrino-free neutron stars. 
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I. INTRODUCTION 

Studying the properties of neutron stars allows us to 
test our knowledge on the properties of nuclear mat- 
ter under extreme conditions. Theoretical studies have 
shown that a neutron star is expected to have a liquid 
core surrounded by an inner crust [l| , which extends out- 
ward to the neutron drip-out region. While the neu- 
tron drip-out density pout has been relatively well deter- 
mined 1^, the transition density pt at the inner edge of 
the crust is still quite uncertain because of our limited 
knowledge on the nuclear equation of state (EOS), es- 
pecially the density dependence of the symmetry energy 
(i?sym(p)) of neutron-rich nuclear matter Recently, 
significant progress has been made in constraining the 
EOS of neutron-rich nuclear matter using terrestrial lab- 
oratory experiments (See Ref. Q for a recent review). 
In particular, from analyses of experimental data on 
neutron skin thickness, isobaric analogue states. Pygmy 
dipole resonances, and giant dipole resonances in nuclei 
as well as on isospin diffusion, isoscaling, and neutron- 
proton to triton-^He ratio in intermediate-energy nuclear 
reactions, significant constraints on Egymip) have been 
obtained for the same sub-saturation density region as 
expected in the inner edge of neutron star crusts. The 
extracted slope parameter L = 3po(9£'sym(p)/5p)p=p(, of 
the nuclear symmetry energy from these studies has val- 
ues in the range 30 MeV < L < 80 MeV [l]. With 
the MDI interaction together with the value L = 86 ± 25 
MeV constrained from an analysis of the isospin-diffusion 
data 0413 in heavy -ion collisions using the isospin- 
dependent Boltzmann-Uehling-Uhlenbeck (IBUU) trans- 
port model with the momentum-dependent MDI interac- 
tion the density and pressure at the inner edge of 
the crust of cold neutron stars were studied by consid- 



ering the boundary of the instability region or the spin- 
odal boundary between the liquid core and inner crust 
of a cold neutron star in both the thermodynamical ap- 
proach 0, [13] and the dynamical approach |13l - [l7| . This 
leads to the constraints 0.040 fm^^ < pt < 0.065 fm~'^ 
and 0.01 MeV/fm^ < Pt < 0.26 MeV/fm^, respectively, 
for the transition density and pressure. Together with 
the crustal fraction of the total moment inertia of the 
Vela pulsar extracted from its glitches [l^ , a tighter con- 
straint on the mass-radius relation of cold neutron stars 
was obtained [l9| . 

Because of the initial high temperature and apprecia- 
ble proton fraction in a newly-formed neutron star imme- 
diately after gravitational collapse of a massive star [20l - 
[22j . neutrinos are abundantly produced from the Urea 
process in its inner core. Although high energy neutrinos 
can be trapped at densities as low as 10^^ g/cm'^ [2^, the 
stars cools by neutrino emissions. As neutrinos are emit- 
ted from this so-called proto-neutron star, which has an 
initial temperature of - lO^^K (about 10 MeV) [Mill, 
its temperature drops to ~ 10^"/^ (about 1 MeV) and 
even lower. Afterwards, the neutron star becomes trans- 
parent to neutrinos as their mean free path increases with 
decreasing energy, and the cooling of the neutron star 
continues to be dominated by neutrino emission for a 
long time. It is thus of interest to study the transition 
density and pressure in newly-born hot neutron stars, as 
this would help to understand the cooling mechanism and 
structural evolution of neutron stars. In this paper, we 
extend the study of Ref. [l^ to finite temperature and 
study the dependence of the transition density and pres- 
sure of hot neutron stars on the nuclear symmetry energy, 
particularly its slope at nuclear saturation density. 

This paper is organized as follows. We first review 
in Sec. In] the momentum-dependent MDI interaction for 
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nucleons, in Sec. IIIII the properties of hot neutron star 
matter, and in Sec. IIVI the dynamical approach for lo- 
cating the inner edge of the crust of a hot neutron star. 
We then show in Sec. |V]the results and conclude with a 
summary in Sec. IVII 



II. THE MDI INTERACTION 

The MDI interaction is an effective nuclear interaction 
with its density and momentum dependence constrained 
from the phenomenological finite-range Gogny interac- 
tion In the mean- field approximation, the poten- 
tial energy density of a nuclear matter of density p and 
isospin asymmetry 5 = {pn Pp)/ P , with p„ and pp be- 
ing, respectively, the neutron and proton densities, can 
be expressed as H, [lH 
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In the above equation, t(t') is the nucleon isospin taken 
to be 1/2 for neutron and —1/2 for proton; /r('^,p) = 
(2//i3){exp[(pV2m + Ur- Pt)/T] + l}-^ is the nucleon 
phase-space distribution function in a thermally equili- 
brated nuclear matter with m — 939 MeV being the nu- 
cleon mass, pr being the chemical potential of nucleon 
of isospin r, T being the temperature and U-r being the 
nucleon mean-field potential to be introduced below; and 
po = 0.16 fm~^ is the saturation density of normal nu- 
clear matter. Values of the parameters A„(x), Ai{x), 
B, a, A, Ci = Cr,T and C„ = Ci-,_r can be found in 
Refs. @, For symmetric nuclear matter, this inter- 
action gives a binding energy of —16 MeV per nucleon 
and an incompressibility Kq of 212 MeV at saturation 
density. 

Taking the derivative of Eq. ([1]) with respect to the 
proton or neutron density leads to the following single- 
particle potential for a nucleon of isospin r: 

U{p,S,p,T) = Au{x)^ + Ai{x)^ 
Po Po 

+ B (^-^y {1 ^ xs') 

— 8tx — (5p-T 

<y + i P^o 

Po J 



Po 



l + (j7-p')2/A2 
,3 , f-rir^f) 



l + (p-p')2/A2 



(2) 



which is seen to depend on the momentum p of the nu- 
cleon. Although the properties of cold nuclear matter can 
be essentially determined analytically using Eqs. ([T]) and 
([5]), these equations needs to be solved numerically and 
self-consistently by the iteration method [2^ to obtain 
the thermodynamical quantities of hot nuclear matter. 

In cold nuclear matter, the symmetry energy from the 
MDI interaction is given by 



Esym{p) 



1 f d^E 

2 [~d5^ 



+ ^(-24.59 + 4B./(. + l)) 



Bx 
a+1 \po 



4pj - A^pj In 



C, 



47r 



9pop \ 
fApj + A^^ 
\ A2 



A^ 



9pop \ J 
4pj-p}i'ip} + A'')ln 



'Apj -I- A^' 
A2 



(3) 



where pj = h{3T:^ p/2)^^^ is the nucleon Fermi momen- 
tum in symmetric nuclear matter. The first term in 
Eq. ([3]) is the contribution from the kinetic part while 
other terms are from the potential part. The symmetry 
energy is fixed to be 30.5 MeV at normal nuclear density, 
and the parameter x is used to model the density depen- 
dence of the symmetry energy away from the saturation 
density without changing the properties of symmetric nu- 
clear matter. The resulting slope parameter of the sym- 
metry energy has values of about 15 MeV, 60 MeV, and 
106 MeV for a; = 1. 0; and —1, respectively. 

The density dependence of the symmetry energy from 
the MDI interaction is shown in panel (a) of Fig. [1] for 
X = I, X ~ and x = —1. With x = 1 {x = —1) the sym- 
metry energy is larger (smaller) at subsaturation densi- 
ties but smaller (larger) at suprasaturation densities, and 
it is called a 'soft' ('stiff') symmetry energy. For a; = 0, 
the value of the symmetry energy lies between those of 
X = I and a; = — 1. In panels (b) and (c), the momentum 
dependence of the nucleon single-particle potential Uq in 
symmetric nuclear matter at half saturation density is 
displayed. Also shown in the figure is the symmetry po- 
tential, defined as Usym = {Un — Up)/25 and calculated at 
S = 0.2 and the same density. It is seen that Uq has nega- 
tive values at low momenta but increases with increasing 
momentum and saturates at high momenta. The sym- 
metry potential J/sym, which measures the difference be- 
tween the neutron and proton single-particle potentials 
in an asymmetric nuclear matter, decreases with increas- 
ing momentum and is larger for a; = 1 than for a; = 
and a; = — 1, which is consistent with the behavior of the 
symmetry energy at subsaturation densities. 
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FIG. 1: (Color online) (a) Density dependence of symmetry 
energies, (b) single-particle potentials in symmetric nuclear 
matter and (c) symmetry potentials in cold nuclear matter at 
po/2 for the MDI interaction. 



III. HOT NEUTRON STAR MATTER 
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FIG. 2: (Color online) The proton fraction x^, as a function of 
baryon density p from the MDI interaction for both neutrino- 
trapped matter (left panel) and neutrino-free matter (right 
panel) with different x parameters and at different temper- 
atures. Note that different scales for used for the 
neutrino-free and the neutrino-trapped matter. 



For a newly-born neutron star, its cooling is dominated 
by the emission of neutrinos, which are produced through 
the Urea process, leading to the following charge neutral- 
ity and /3-stable conditions: 
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The production of muons is neghgible for the density and 
temperature present in the hot neutron star crust. Neu- 
trino trapping in the early stage of a supernova has been 
extensively studied in the literature (23| , and it was found 
that the fraction of leptons 
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is about 0.35 ^ 0.4 at the onset of trapping [2l|, |2J| 
and decreases as neutrinos leave the star. In the present 
study, we consider the two extreme cases of neutrino- 
trapped and neutrino-free hot neutron star matters. In 
the former case, we choose Yj to be 0.4 as an example and 
consider the typical temperature of 5 or 10 MeV. In the 
case of neutrino-frec hot neutron star matter, we set the 
neutrino chemical potential to be and the temperature 
to be or 1 MeV, corresponding to the later stage of a 
neutron star's evolution. The abundance of each species 
in both cases can be calculated from above equations to- 
gether with the baryon number conservation condition 

P = P71+ Pp- 

Assuming that both electrons and neutrinos arc mass- 
less, we determine the proton fraction Xp = {\ — S)/2 in 



both the neutrino-trapped matter and the neutrino-free 
matter as a function of baryon density p, and the results 
are shown in Fig. [2] It is seen that the neutrino- free mat- 
ter is much more neutron-rich than the neutrino-trapped 
matter. The critical proton fraction 11 ^ 15 % for the 
direct Urea process [13, [l^ is smaller than the proton 
fraction in the neutrino-trapped matter but larger than 
the proton fraction in the neutrino-frec matter. The pro- 
ton fraction increases with increasing density in all cases 
and slightly increases with increasing temperature for a 
fixed density, especially at low densities. These can be 
understood from the increase of the symmetry free en- 
ergy with increasing temperature and density [26| , which 
makes the neutron star matter more symmetric at higher 
temperatures and densities. Furthermore, the stiff sym- 
metry energy (x = —1) makes the system more neutron- 
rich at subsaturation densities as expected. 

The total pressure P in a hot neutron star matter can 
be written as 



(7) 



where the baryon contribution Pb{p,T,S) is calculated 
from the thermodynamic relation 



Pbip,T,d) = 



TY,Sr{p,T,S)~Vip,T,S) 

T 
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In the above equation, V{p,T,S) and 14;in(p, T, (5) are, 
respectively, the potential and kinetic contributions to 
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the total energy density with the latter given by 

Vkin [p, T,S)^J2 [ d^P^fr (r, P), (9) 
^ — ' / Zm 

T 

and St-(p, T, S) is the entropy density, which is given by 

Sr{p,T,6) = -j-^ / p'^[nrlnnr + il-nr)\n{l-nr)]dp, 

(10) 

with the particle occupation number 

(11) 



exp[(p2/2m + Ur- fir)/T] + 1 ■ 



Above formula can also be used for the calculation of P/ 
by using I (l =■ e, z^e) instead of t, and leptons are treated 
as non-interacting ultra-relativistic particles. 



IV. LOCATING THE INNER EDGE OF THE 
NEUTRON STAR CRUST 

We briefly review in this section the dynamical ap- 
proach that will be used to locate the inner edge of the 
neutron star crust and discuss its application to the 
case of finite temperature. We neglect muons in the neu- 
tron star crust as their number is small compared to that 
of electrons at low densities and temperatures. 

In the dynamical approach, the stability condition for a 
homogeneous neutron star matter against small periodic 
density perturbations can be well approximated by [isl - 

[a 



> 0, (12) 



where k is the wavevector of the spatially periodic density 
perturbations and 
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The three terms in Eq. (|12l) represent, respectively, con- 
tributions from the bulk nuclear matter, the density- 
gradient terms, and the Coulomb interaction. The empir- 
ical values for the coefficients of density-gradient terms 
are Dpp = L>„„ = Dnp = 132 MeV-fm^ [H [H. At 
fcmin = [(47re2//3)^/2 — k'^pY^'^., Vdyn(fc) has the minimal 
value of ydyn(fcmin) ^ + 2{A^e'Pf'^ - Pklp dMl, 
and pt is then determined from Vdyn(fcmin) = 0. We note 
that the first term in Eq. ([12]) gives the dominant con- 
tribution in the determination of the transition density, 
and including other terms lowers the transition density. 
Also, although at low temperatures the electron density 



Pe can be written as an expansion of the electron chemi- 
cal potential 
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which gives an analytical expression for dpe/dpe, at high 
temperatures numerical calculations are needed. 

The dynamical approach reduces to the so-called ther- 
modynamical approach [J, [l^ in the long wave length 
limit when the density gradient terms and the Coulomb 
interaction are neglected [H, [ll] , which leads to the sta- 
bility condition 
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FIG. 3: (Color online) The instability region and the rela- 
tive neutron-proton abundance of hot neutron star matter for 
different temperatures and nuclear symmetry energy param- 
eters. Regions where they cross each other are shown in the 
insets in enlarged scales. 

To illustrate the relation between the transition density 
and the area of the spinodal region, we show in Fig.[3]the 
instability region of nuclear matter with the boundary 
determined by Vthcr = and the relative neutron-proton 
abundance of a hot neutron star in the {pmPp) plane. 
The cross point, also shown with enlarged scales in the 
insets, is the transition density from the thermodynam- 
ical approach. It is seen that although the neutron star 
matter becomes less neutron-rich with increasing tem- 
perature, the area of the instability region shrinks more 
quickly with increasing temperature. Furthermore, the 
stiffness of the symmetry energy also affects the shape 
and area of the spinodal region. As temperature in- 
creases, the spinodal boundaries from different values of x 
cross with the relative neutron-proton abundance curves 
at decreasingly small proton densities. For a more de- 
tailed discussion on the symmetry energy and tempera- 
ture effects on the spinodal region, we refer readers to 
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Refs. [30H35|- The following analysis of the transition 
density and pressure in hot neutron stars is, however, car- 
ried out by using the more realistic dynamical approach. 



V. RESULTS AND DISCUSSIONS 

In this section, we show the temperature dependence 
of the transition density and pressure in newly-born hot 
neutron stars by using the MDI interaction with different 
values for the symmetry energy parameter x or the slope 
parameter L of the symmetry energy. 
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FIG. 4: (Color online) Transition densities pt ((a) and (c)) 
and pressure Pt ((b) and (d)) as functions of the slope pa- 
rameter L of the symmetry energy at different temperatures 
for both the neutrino-trapped matter ((a) and (b)) and the 
neutrino- free matter ((c) and (d)). 

The dependence of the transition density and pressure 
in hot neutron stars on the slope parameter L of the 
symmetry energy at different temperatures is shown in 
Fig. m It is seen that the transition density pt generally 
decreases with increasing value of the slope parameter 
L of the symmetry energy. As the transition density can 
be viewed approximately as the beginning of a first-order 
liquid-gas phase transition, a stiffer symmetry energy, 
which corresponds to a softer equation of state at subsat- 
uration densities, leads thus to a smaller phase transition 
density and therefore a lower core-crust transition den- 
sity. Furthermore, the transition density decreases with 
increasing temperature, and for the neutrino-free matter 
this is more pronounced for larger values of L. The tem- 
perature effect can again be understood by the decreas- 
ing phase transition density with increasing temperature. 
For the neutrino-trapped matter, the L-dependence of 
the transition density is relatively weak. The weak L- 
dependence is mainly due to the fact that the isospin 
asymmetry is not so large as shown in Fig. [5J We note 
that a similar temperature dependence of the transition 
density has been obtained in studies based on Skyrme in- 



teractions and relativistic mean field models [3J, [361, l37| . 
For the L-dependence of the transition density, results 
from these models are, however, not so clear as differ- 
ent values for the incompressibility Kq of the symmetric 
matter at saturation density and Esymipo)^ which also 
affect the value of the transition density, have been used. 

For the transition pressure Pt , its value in the neutrino- 
trapped matter decreases only very slightly with increas- 
ing L as a result of the weak L-dependence of pt- Its 
temperature dependence shows, however, a complicated 
behavior of slightly higher and smaller values at higher 
temperatures for smaller and larger values of L, respec- 
tively. This is due to the fact that although the tran- 
sition density pt decreases with increasing temperature, 
the contribution from leptons increases with increasing 
temperature. For the neutrino-free matter, Pt is seen to 
decrease rapidly with increasing L. As to its temperature 
dependence, Pt in the neutrino-free matter decreases with 
increasing temperature at larger values of L but shows a 
weaker temperature dependence for smaller values of L. 
Also, Pt is larger for the neutrino-trapped matter than 
for the neutrino-free matter. Interestingly, Pt becomes 
very small and even negative at T = f with larger L for 
the neutrino-free matter. This is due to the smaller con- 
tributions from the leptons and the asymmetric part of 
the nuclear interactions to the total pressure. Since the 
pressure at the inner edge of neutron star crust cannot 
be negative, our finding thus indicates that either the 
neutrino-free matter in hot neutron stars cannot reach a 
temperature above T = 1 MeV or the symmetry energy 
cannot have a slope parameter larger than L ^ ICQ MeV. 
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FIG. 5; (Color online) Transition densities pt ((a) and (c)) 
and pressure Pt ((b) and (d)) as functions of temperature 
T with X = and x = —1 for both the neutrino-trapped 
matter ((a) and (b)) and the neutrino- free matter ((c) and 
(d)). Note that different scales for temperatures are used for 
the neutrino-trapped matter and the neutrino-free matter. 

The temperature effect on the transition density and 
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pressure for a fixed symmetry energy parameter is 
demonstrated in Fig. [5] for the symmetry energy param- 
eters X = 1, X — and x ~ —1. For the neutrino- 
trapped matter, the temperature effect is similar for all 
three x values, and the transition density pt decreases al- 
most linearly with increasing temperature at lower T and 
decreases quickly at higher T. For the neutrino-free mat- 
ter, although the transition density pt decreases smoothly 
with increasing temperature for a; = 1 , it decreases slowly 
(quickly) at lower temperatures with a; = (a; = — 1) 
while quickly (slowly) at higher temperatures. The tem- 
perature effects on the transition density reflect those on 
the spinodal region and the abundance of particle species 
in the hot neutron star matter. For the neutrino-trapped 
matter, the transition pressure Pt is seen to be insensitive 
to the temperature for x = — 1 but increases slightly with 
increasing temperature for x ^ and .t = 1 at lower T, 
and it decreases with increasing temperature at higher 
T for all values of x. For the neutrino-free matter, it is 
insensitive to temperature for x = 1 but decreases with 
increasing temperature for a; = 0, while for a; = — 1 it 
drops to a negative value and becomes positive again as 
the temperature increases. Our results again show that 
the behavior of Pt is dominated by that of pt and the 
contribution from the leptons, with the former decreasing 
and the latter increasing with increasing temperature. 
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temperature Tc at which the transition density pt dis- 
appears is thus useful for understanding the structural 
evolution of newly-born hot neutron stars. Figure [5] dis- 
plays the L-dependence of the critical temperature for 
both the neutrino-trapped matter and the neutrino-free 
matter. One sees that the critical temperature Tc de- 
creases slightly with increasing value of L in the neutrino- 
trapped matter, but it first decreases and then increases 
with increasing L in the neutrino-free matter. This com- 
plicated behavior is due to the isospin and temperature 
effects on the spinodal region and the relative neutron- 
proton abundance as shown in Fig. [31 Our results thus 
indicate that for neutrino-trapped neutron stars of tem- 
peratures higher than 12 MeV or for neutrino-free neu- 
tron stars of temperatures higher than 1.5 MeV, there 
exists no inner crust if the value of L is 70 ^ 80 McV. 
As newly-born hot neutron stars cool, the temperature 
at which the inner crust can form thus depends on the 
density dependence of the symmetry energy at subsat- 
uration densities. Again, the magnitude of the critical 
temperature for both neutrino-trapped and neutrino-free 
matter is similar to those from Skyrme interactions and 
relativistic mean field models [13, 1131 • 

The above results were obtained using the dynamical 
approach that includes the effects of both the density 
gradient terms and the Coulomb interaction. Neglecting 
these effects, the resulting thcrmodynamical approach 
given by Eq. (1151) gives higher values for both the transi- 
tion density and pressure. This is especially the case for 
the neutrino-trapped matter and/or for smaller values of 
L as more electrons are present in such hot neutron stars 
and the effect due to the Coulomb interaction becomes 
more important. Also, there are recently some studies 
on the transition density in neutron stars using various 
nucleon-nucleon interactions [H, H^. To see the effect 
of momentum dependence in the nucleon-nucleon inter- 
action on the transition density and pressure, we have 
also repeated above calculations using the momentum- 
independent MID interaction [4^, which gives the same 
equation of state for asymmetric nuclear matter but dif- 
ferent single-particle mean-field potential in comparison 
with the momentum-dependent MDI interaction used in 
the present study, and we find that the momentum- 
dependent effect on the transition density and pressure 
in hot neutron stars is small. 



FIG. 6: (Color online) Critical temperature Tc as a function 
of the slope parameters L of the symmetry energy for both 
the neutrino-trapped matter and the neutrino-free matter. 

We have seen in Fig. [3] that as the temperature of the 
neutron star matter increases, there will be eventually 
no cross point between the curve of neutron-proton rela- 
tive abundance and the boundary of the spinodal region, 
leading to the disappearance of the transition density in 
hot neutron stars. In such case, the inner crust (nuclear 
'pasta' phase) disappears and the liquid core expends 
directly to the outer crust. To determine the critical 



VI. SUMMARY 

We have studied the transition density and pressure 
at the boundary that separates the liquid core from 
the inner crust of neutron stars using the momentum- 
dependent MDI interaction in both the neutrino-trapped 
matter and the neutrino-free matter at finite tempera- 
tures, which are expected to exist during the early evo- 
lution of neutron stars. In particular, we have investi- 
gated the effect of nuclear symmetry energy by varying 
the parameter x in the MDI interaction from 1 to — 1, 
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corresponding to the values 15 < L < 106 MeV for the 
slope of nuclear symmetry energy at normal density that 
were constrained by both the isospin diffusion data 0-01 
and other experimental observables @. We have found 
that the transition density and pressure are larger in the 
neutrino-trapped matter than in the neutrino-free mat- 
ter. Furthermore, the transition density and pressure are 
found to roughly decrease with increasing temperature 
and L for both the neutrino-trapped and the neutrino- 
free matter, except that the transition pressure shows a 
complicated relation to the temperature for the neutrino- 
trapped matter. Also, negative values of the pressure at 
the transition density have been obtained, which can be 
used to rule out a very stiff symmetry energy at sub- 
saturation densities. We have also studied the critical 
temperature above which the inner crust (nuclear 'pasta' 
phase) cannot be formed in newly-born neutron stars and 
found that it depends sensitively on the density depen- 
dence of the nuclear symmetry energy at subsaturation 



densities. 
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